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ABSTRACT 

c/^ : 

! We use joint likelihood analyses of combinations of fifteen cosmic microwave back- 

; ground (CMB) anisotropy data sets from the DMR, UCSB South Pole 1994, Python 

^ ; I-III, ARGO, MAX 4 and 5, White Dish, OVRO, and SuZIE experiments to constrain 

' cosmogonies. We consider open and spatially-flat-A cold dark matter cosmogonies, 

■ with nonrelativistic-mass density parameter in the range 0.1-1, baryonic-mass den- 
Q\ • sity parameter fi^ in the range (0. 005-0. 029)/i~^, and age of the universe to in the range 
O : (10-20) Gyr. 

, Marginalizing over all parameters but the data favor 0,q ~ 0.9-1 (0.4-0.6) flat- 

r~| ! A (open) models. The range in deduced Qq values is partially a consequence of the 

different combinations of smaller-angular-scale CMB anisotropy data sets used in the 
O ' analyses, but more significantly a consequence of whether the DMR quadrupole moment 

' is accounted for or ignored in the analysis. While the open model is difficult to reconcile 

■ with the results of less exact analyses of more recent CMB anisotropy data, the lower 
J> . values of found in this case are more easily reconciled with dynamical estimates 

^ , of this parameter. For both flat-A and open models, after marginalizing over all other 

^ [ parameters, a lower Q.sh'^ — 0.005-0.009 is favored. This is also marginally at odds with 

estimates from more recent CMB anisotropy data and some estimates from standard 
nucleosynthesis theory and observed light element abundances. For both sets of models 
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a younger universe with to — 12-15 Gyr is favored, consistent with other recent non- 
CMB indicators. We emphasize that since we consider only a small number of data 
sets, these results are tentative. More importantly, the analyses here do not rule out 
the currently favored flat-A model with ~ 0.3, nor the larger QbIt-^ values favored 
by some other data. 

Subject headings: cosmic microwave background — cosmology: observations — large-scale 
structure of the universe 



1. Introduction 

There has been a remarkable increase in the quality and quantity of cosmic microwave back- 
ground (CMB) anisotropy measurements since the initial detection of the anisotropy on large an- 
gular scales a decade ago.^ These measurements are becoming increasingly useful in the continuing 
processes of determining how well cosmological models approximate reality and for constraining 
cosmological parameters such as JIq, h, and CIb in these models^. 

In this paper we utilize the full information in the data from each experiment in an effort 
to place robust constraints on cosmological parameters. This is achieved through a maximum 
likelihood analysis of the data using realistic model anisotropy spectra. Most of the data sets we 
consider in this paper are small, in contrast to some of the more recent data sets (see, e.g., Netterfield 
et al. 2002; Stompor et al. 2001), which because of their size require a more approximate analysis 
technique. Ganga et al. (1997a, hereafter GRGS) extend the maximum likelihood technique to 
account for uncertainties, such as those in the beamwidth of the telescope and the calibration of 
the experiment.^ This technique has been used with model CMB anisotropy spectra in analyses of 
the Gundersen et al. (1995) UCSB South Pole 1994 data, the Church et al. (1997) SuZIE data, 
the Lim et al. (1996) MAX 4+5 data, the Tucker et al. (1993) White Dish data, the dc Bcrnardis 
et al. (1994) ARGO data, the Piatt et al. (1997) Python I-III data, and the Leitch et al. (2000) 
OVRO data (GRGS; Ganga et al. 1997b, 1998; Ratra et al. 1998, 1999a, hereafter R99a; Rocha et 
al. 1999; Mukherjee et al. 2002b). 

Given the error bars associated with these measurements, interesting constraints on cosmo- 



''See, e.g., Miller et al. (2002a), Coble et al. (2001), Scott et al. (2002), and Mason et al. (2002) for recent 
measurements. 

*Here h is the Hubble constant -ffo in units of 100 km s^^ Mpc^^. See, e.g., Podariu ct al. (2001), Wang, Tegmark, 
& Zaldarriaga (2002), Miller et al. (2002b), Sievers et al. (2002), and Lewis & Bridle (2002) for recent discussions of 
constraints on cosmological parameters. 

®For some of the data sets we consider in this paper foreground non-CMB contamination must also be accounted 
for; see, e.g., Kogut et aJ. (1996), de Oliveira-Costa et al. (1998), Hamilton & Ganga (2001), and Mukherjee et al. 
(2002a, 2002b) for discussions of the method used to accomplish this. 
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logical model parameters require the joint analysis of many data sets. If the measurements are 
acquired for regions well separated in space, or on very different angular scales, the likelihoods 
of the individual data sets arc independent and can thus be multiplied together to construct the 
likelihood of the combined data set. This combined likelihood is then used to derive constraints on 
cosmological model parameters. A combined analysis of the smaller-angular-scale data sets listed 
above, excluding the Python and OVRO data, is presented in Ratra et al. (1999b, hereafter R99b). 
In this paper we extend the analysis of R99b to include the Python and OVRO data, as well as the 
large angular scale DMR data (Gorski et al. 1998, hereafter G98; Stompor 1997). 

In § 2 we describe the models and cosmological parameter space we consider. See R99a for a 
more detailed description. In § 3 wc summarize the various combinations of data sets we consider. 
See R99b for further details. In § 4 we summarize the computational techniques we use. See GRGS 
and R99b for more detailed discussions. In § 5 we present and discuss results from analyses of 
various combinations of the smaller-angular-scale CMB anisotropy data sets, and in § 6 we add the 
DMR data to the mix. We conclude in § 7. 

2. Cosmogonical Models 

In this paper we focus on a spatially-flat cold dark matter (CDM) model with a cosmological 
constant A. This model is consistent with most current measurements.^'^ While it is important to 

compare the currently favored model to the CMB anisotropy data, it is also important to check 
how well other, currently less favored, models fare. We therefore also examine the constraints these 
data place on a spatially-open model with no A (see, e.g., Gott 1982; Ratra k. Peebles 1995). For 
more detailed discussions of these models see R99a and R99b. 

The CMB anisotropy spectra in these models are generated as quantum fluctuations in a 
weakly coupled field during an early epoch of inflation. They are thus realizations of spatially 
stationary Gaussian random processes (see, e.g., Ratra 1985; Fischler, Ratra, & Susskind 1985). 
The measured CMB anisotropy appears to be consistent with this Gaussianity assumption (see, 
e.g., Mukherjee, Hobson, & Lasenby 2000; Park et al. 2001; Wu et al. 2001; Shandarin et al. 2002; 
Santos et al. 2002; Komatsu et al. 2002; Polenta et al. 2002), and the experimental noise also 
appears to be Gaussian, thus validating our use of the likelihood analysis method of GRGS. 

To make the computations tractable, in each model (flat-A and open) we consider CMB 
anisotropy spectra parameterized by (i) the quadrupole-moment amplitude Qrms-PS; (ii) f^Oj (iii) 



^°See, e.g., Peebles (1984), Efstathiou, Sutherland, & Maddox (1990), Stompor, Gorski, & Banday (1995), Ratra 
et al. (1997), Sahni & Starobinsky (2000), Carroll (2001), and Peebles & Ratra (2002) for discussions of this model. 
While not considered in this paper, a spatially-flat model dominated at the current epoch by time-variable dark 
energy is also largely consistent with current measurements (see, e.g., Peebles & Ratra 1988; Ratra & Quillen 1992; 
Steinhardt 1999; Brax, Martin, & Riazuelo 2000; Munshi & Wang 2002; Bartelmann, Perrotta, & Baccigalupi 2002; 
Dave, Caldwell, & Steinhardt 2002; Chen & Ratra 2002; Podariu et al. 2002). 
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^Bh^, and (iv) tg. While it is of interest to also consider other cosmological parameters current 
data do not require consideration of a larger dimensional parameter space. That is, current data 
appear to indicate that effects like reionization and tilt and those due to primordial gravity waves 
are small. Furthermore, the maximum likelihood analysis technique used here is computationally 
intensive and an analysis that takes account of more cosmological parameters will be very time 
consuming. Rather than make approximations (e.g., data compression) to speed up such a compu- 
tation and minimize memory requirements, in this paper we perform as accurate a computation as 
possible based on model spectra from an as small as viable yet still realistic cosmological parameter 
space. 

The CMB anisotropy spectra we use are computed for a range of $7o spanning the interval 0.1 
to 1 in steps of 0.1, for a range of Qsh"^ spanning the interval 0.005 to 0.029 in steps of 0.004, and 
for a range of io spanning the interval 10 to 20 Gyr in steps of 2 Gyr. In total 798 spectra are 
computed to cover the cosmological-parameter spaces of the open and fiat-A models. Examples of 
spectra are shown in Fig. 2 of R99a, Fig. 1 of R99b, and Fig. 2 of Rocha et al. (1999). 

3. CMB Anisotropy Data Sets 

R99b consider various combinations of ten different data sets. These were the UCSB South Pole 
1994 Ka and Q band observations, hereafter SP94Ka and SP94Q (Gundersen et al. 1995; GRGS), 
the ARGO Hercules observations (de Bernardis et al. 1994; R99a), the MAX 4 i Draconis (ID) 
and o Hercuhs (SH) and MAX 5 HR5127 (HR), /i Pegasi (MP), and (j) Herculis (PH) observations 
(Lim et al. 1996; Ganga et al. 1998), the White Dish observations (Tucker et al. 1993; Ratra et al. 
1998), and the SuZIE observations (Church et al. 1997; Ganga et al. 1997b). Detailed information 
about these data sets may be found in the papers cited above and in R99b. 

In this paper we augment the data sets used in R99b with the Python I-III data (Piatt et 
al. 1997; Rocha et al. 1999), the OVRO data (Leitch et al. 2000; Mukherjee et al. 2002b), and 
the DMR data (Bennett et al. 1996; G98; Stompor 1997). This adds an additional five CMB 
anisotropy data sets to the ten considered in R99b. 

R99b consider three different combinations of the CMB anisotropy data sets. The first com- 
bination included all data (SP94, ARGO, MAX 4 and 5, White Dish, and SuZIE). The second 
combination excluded the SuZIE data, which probe the smallest angular scales (to multipolc index 
I ~ 2000 in some models, Ganga et al. 1997b), since our model CMB anisotropy spectra do not 
account for effects that could be important on these small angular scales. The third combination of 
data used in R99b include just the measurements of SP94Ka, MAX 4 ID, and MAX 5 HR which are 
thought to provide the most reliable constraints, among the data considered in R99b (see discussion 
in R99b) , to constrain the parameters of the theoretical CMB anisotropy spectra considered in that 
paper. 

The new smaller-angular-scale data used in this paper, Python I-III and OVRO, raise no 
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obvious red flags : we have no reason to think that they are any less rehable than the data 

considered to be the most rehable in R99b. We therefore again derive cosmological constraints for 
three different combination CMB anisotropy data sets, augmenting each of the three combination 
data sets of R99b with the Python I-III and OVRO data. 

The situation with the largest angular scale CMB anisotropy data, that measured by the 
DMR experiment (Bennett et al. 1996; G98), is more complicated. While a couple of other 
issues also need to be considered, the main complication arises from the effect on the derived 
cosmological constraints from the inclusion or exclusion of the observed DMR quadrupole moment 
(Kogut et al. 1996) in the analysis (G98; Stompor 1997). Non-CMB anisotropy emission from the 
Milky Way is predominantly quadrupolar, and any post-correction remnant of this will primarily 
affect the observed DMR quadrupole moment. To be sure, this is only of order a 1 cr effect on 
cosmological parameter determination from the DMR data alone. However, given the differences 
in the constraints derived from the different combinations of smaller-angular-scale data sets we 
consider (and again, these differences are of order 1 cr), folding in the DMR data results in more 
significant differences between cosmological constraints derived using the various combination data 
sets. One might hope that better maps of foreground emission will help resolve the DMR quadrupole 
moment issue, but that is for the future. The only option at present appears to be to consider 
cosmological constraints that follow on using all "reasonable" combinations of data. To illustrate 
the situation we consider the two most extreme DMR data combinations (G98): (i) the galactic 
frame data including the quadrupole moment and correcting for faint high-latitude galactic emission 
(this results in the smallest Qrms-PS value); and (ii) the ecliptic frame data with the quadrupole 
moment excluded and no other galactic emission correct ion (this results in the largest (^rms— PS 
value). In conjunction with the three different combinations of smallcr-angular-scale data, these 
two DMR data sets result in six different combinations of CMB anisotropy data we consider in this 
paper. 



4. Summary of Computation 

ORGS describe the computation of the likelihood function for a given CMB anisotropy data 
set. Offsets and gradients removed from the data are accounted for in the analysis. Beamwidth 
and calibration uncertainties are also accounted for as described in GRGS. In some cases non-CMB 
anisotropy foreground contamination have also been accounted for. 

The data sets considered in this analysis were acquired from regions well separated in space, 
or are at very different angular resolution. Consequently they are statistically independent and 
the likelihood functions of the individual data sets are simply multiplied together to construct the 



^^We note that Mukherjee et al. (2002b) account for the removal of non-CMB foreground contamination in the 
OVRO analysis. 
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likelihood function of the combined data. The fiat-A and open model hkehhoods we consider are 

a function of the four parameters described in § 2: Qrms-PS! ^O) ^sh'^, and to- MarginaHzed 
hkehhood functions are derived by integrating over one or more of these parameters. We assume a 
uniform prior in the parameters integrated over, set to zero outside the parameter range considered. 

To determine central values and limits from the likelihood functions we assume a uniform 
prior in the relevant parameter. Then the corresponding posterior probability density distribution 
function vanishes outside the chosen parameter range and is equal to the likelihood function inside 
this range. The deduced central value of the parameter is taken to be the value at which the 
posterior probability density peaks, and we quote highest posterior density limits. See GRGS, 
R99a, and R99b for details. The quoted limits depend on the prior range considered for the 
parameter. This is a significant effect if the likelihood function is not sharply peaked within the 
parameter range considered, as is the case for a number of the likelihood functions derived in this 
paper. See R99a and R99b for more detailed discussions of this issue. In what follows we consider 
1, 2, and 3 a highest posterior density limits which include 68.3, 95.4, and 99.7% of the probability. 

Rocha et al. (1999) and Mukherjee et al. (2002b) note that the four-dimensional posterior 
probability density distribution function L(Q]-ms-PS) ^O) ^s^^i ^o) is nicely peaked in the Qrms-PS 
direction but fairly flat in the other three directions. Marginalizing over Qrms-PS results in a 
three-dimensional posterior distribution L(r2o, to) which is steeper, but still relatively flat. 

As a consequence, limits derived from the four- and three-dimensional posterior distributions are 
generally not highly statistically significant. We therefore do not show contour plots of these 
functions here. Marginalizing over Qrms-PS and one other parameter results in two-dimensional 
posterior probability distributions which are more peaked. See Fig. 1 for examples of the constraints 
from the two-dimensional posterior probability distribution functions. As in the ARGO (R99a), 
Python I-III (Rocha et al. 1999), OVRO (Mukherjee et al. 2002b), and earlier combination (R99b) 
data set analyses, in some cases these peaks are at an edge of the parameter range considered. 



5. Smaller Angular Scale Data: Results and Discussion 

Figure 1 shows that even just using the smaller angular scale CMB anisotropy data, the two- 
dimensional posterior distributions allow one to distinguish between different regions of parameter 
space at a fairly high formal level of confidence. For instance, for the SP94Ka, Python I-III, 
MAX 4 ID, MAX 5 HR, and OVRO data combination (panels c and / in the bottom row of 
Fig. 1), the flat-A model near JIq ~ 0.1, JIb/i^ ~ 0.03, and to ~ 10 Gyr, is formally ruled out at 
~ 3 fj confidence. However, we emphasize, as discussed in Mukherjee et al. (2002b), care must 
be exercised when interpreting the discriminative power of these formal limits, since they depend 



^^See Fig. 4 of R99a, Fig. 2 of R99b, Fig. 3 of Rocha et al. (1999), and Fig. 1 of Mukherjee et al. (2002b) for 
related cosmological constraints from subsets of the data considered here. 
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sensitively on the fact that the uniform prior has been set to zero outside the range of the parameter 
space we have considered. 

Figure 2 shows the contours of the two-dimensional posterior distribution for Qrms-PS S'ld CIq, 
derived by marginahzing the four-dimensional distribution over Jl^^^ and t^. These are shown for 
the different combinations of smaller-angular-scalc data and for the DMR data, for both the open 
and flat-A models. Constraints on cosmological parameters from the smaller-angular-scale data are 
consistent with those from the DMR data.^^ It is interesting that the constraints on Qrms-PS derived 
from the DMR data are not very much tighter than those derived from the different combinations 
of smaller-angular-scale data. 

Figure 3 shows the one-dimensional marginalized posterior distribution functions for Clo,Q,Bh^, 
and to- For the flat-A model (first column of Fig. 3), all three combinations of data sets favor JIq = 1, 

and low values of Qq ~ 0.1 are disfavored, at least by those data set combinations that include 
most of the data. This should be compared to the earlier analysis of Fig. 5 of R99b (which did 
not include the Python I-III and OVRO data considered here), where the data were not able to 
significantly discriminate between different fio values in the flat-A case. For the open model (second 
column of Fig. 3) all three combinations of smaller-angular-scale data sets favor 0.4, with low 
values of J7o ~ 0.1 again disfavored. This should also be compared to the results shown in Fig. 5 of 
R99b, where it was found that conclusions about the favored value of f^o in the open model were 
very dependent on the combination of smaller-angular-scale data used to derive them. Apparently 
we now have a large enough collection of data sets for this to no longer be an issue. 

The central two columns of Fig. 3 show the one-dimensional marginalized posterior distribution 
functions for Qsh^- Independent of the data set combination considered low values of Os/i^ are 
favored, with Ush'^ — 0.005 (0.007-0.008) favored in the open (flat-A) model. These results are 
statistically much more significant than, but quite consistent with, those shown in Fig. 6 of R99b. 

The last two columns of Fig. 3 show the one-dimensional marginalized posterior distribution 
functions for to. Independent of the data set combination considered the flat-A case favors an old 
universe with to — 18-20 Gyr (fifth column of Fig. 3), while the open model favors a young universe 
with to ~ 10-12 Gyr (sixth column of Fig. 3). These open model results are very consistent with 
those shown in Fig. 7 of R99b, while R99b found that flat-A models favored a young universe with 
to — 10-12 Gyr, at odds with what we find here. 

More precisely, the one-dimensional distributions for the "most reliable" data in Fig. 3 (the 
bottom row) indicates that an open (flat-A) model with f^o = 0.42 (1.0), or r^^/i^ = 0.005 (0.008), 
or to = 10 (18) Gyr is favored, among the models considered. At 2 a confidence the data formally 
rule out only small regions of parameter space: the data require > 0.18 (Oq > 0.10), or ^Bh^ 
< 0.027 {nBh? < 0.026), or to < 19 Gyr {to > 11 Gyr) for the open (flat-A) model at 2 a. 



"See Fig. 5 of R99a, Fig. 3 of R99b, Fig. 4 of Rocha ot al. (1999) and Fig. 2 of Mukherjee et al. (2002b) for related 
constraints from subsets of the data considered in this paper. 
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While the statistical significance of the constraints on cosmological parameters from these 
combinations of smaller-angular-scale CMB anisotropy data is not high, there arc a couple of 
puzzling issues. The lower Qq ~ 0.4 favored in the open case is more easily reconciled with most 
other estimates of than is the higher value of = 1 favored by the CMB anisotropy data in 
the flat-A model (see, e.g., Peebles Sz Ratra 2002). We emphasize however that more recent CMB 
anisotropy data favor the flat-A model with low JIq (see, e.g., Netterfield et al. 2002; Pryke et 
al. 2002; Stompor et al. 2001; Scott et al. 2002; Mason ct al. 2002). The low values of fis/i^ 
favored by the data considered here are also somewhat at odds with the higher values favored by 
more recent CMB anisotropy data (see, e.g., Netterfield et al. 2002; Pryke et al. 2002; Stompor 
et al. 2001) and those favored by standard nucleosynthesis theory and the observed deuterium 
abundances (Buries, Nollett, & Turner 2001); they are, however, more consistent with the lower 
Q,Bh^ value from standard nucleosynthesis theory and the observed helium and lithium abundances 
(Cyburt, Fields, & Olive 2001). The ages of an open and a flat-A universe determined from the 
smaller-angular-scale CMB data used in this paper lie at the opposite ends of the 2 a range of 
estimates based on the ages of the oldest stars: 11 Gyr ^ to ^ 17 Gyr (sec, e.g., Carrctta ct al. 
2000; Krauss & Chaboyer 2001). Given that the CMB anisotropy data considered here arc not the 
most recent, and that the statistical significance of the puzzling aspects of the derived results are 
not high, perhaps similar analyses of more recent and future data will show that these puzzles are 
only apparent. 

Figure 4 shows the one-dimensional posterior distribution functions for Qrms-PS) derived by 
marginalizing the four-dimensional ones over the other three parameters. The three different com- 
binations of smallcr-angular-scalc data sets considered here result in fairly tight constraints on 
<5rms-PS- At 2 (7 confidence they are consistent with the DMR results for both the open and flat-A 
models. 

The peak values of the one-dimensional posterior distributions shown in Figs. 3 and 4 are 
listed in the figure captions for the case when the four-dimensional posterior distributions are 
normalized such that L((5rms-PS = /xK) = 1. With this normalization, marginalizing over 
the remaining parameter the fully marginalized posterior distributions are, for the open (flat-A) 
models: 1 x 10^^2(4 ^ io27i) for all the data; 1 x 10'^'^'^{2 x lO^^^^) for ah the data cxchiding SuZIE; 
and, 1 x 10^°^ (8 ^ 10^"!) for the SP94Ka, Python I-III, MAX 4 ID, MAX 5 HR, and OVRO 
combination. These numerical values are consistent with the indications from the first two columns 
of Fig. 3 that the data mildly favor the open model over the fiat-A one (since the rig = 1 models 
at the right edges of pairs of panels in each row of Fig. 3 are identical) . 

6. Including the DMR Data: Results and Discussion 

As mentioned above, cosmological parameter constraints derived from the DMR data are some- 
what sensitive to whether the measured DMR quadrupole moment is included in or excluded from 
the analysis. To account for this source of uncertainty we consider the two most extreme (in 
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terms of Qrms-PS normalization) DMR data combinations (G98; Stompor 1997): (i) the galactic 
frame data including the quadrupole moment and correcting for faint high-latitude galactic emis- 
sion; and (ii) the ecliptic frame data with the quadrupole moment excluded and no other galactic 
emission correction. Since we consider three different combinations of smaller-angular-scale data, 
combining these two DMR data sets with the smaller-angular-scale data results in the six different 
combinations of CMB anisotropy data we consider in this section. 

When the smaller-angular-scale data of § 5 are augmented with the DMR data, the four- 
dimensional posterior probability density distribution function L(Qrms-PS) ^Oj ^s^^j ^o) becomes 
much more peaked in the to direction (in addition to being peaked in the Qrms-PS direction) but 
is still somewhat flat in the other two directions. As discussed in § 5, marginalizing over Qrms-PS 
and one other parameter results in two-dimensional posterior probability functions which are more 
peaked. Examples are shown in Fig. 5. Comparing to the corresponding plots (Fig. 1) based on just 
the smaller-angular-scale data, one sees that the inclusion of the DMR data significantly tightens 
the constraints on cosmological parameters. 

For example, Fig. 5 shows that, independent of data set combination considered, flat-A models 
with parameter values near J7o ~ 0.1, ~ 0.03, and to ~ 20 Gyr, and open models with 

flo ^ 0.2, are both formally ruled out at ~ 3 cr confidence. 

Figure 6 shows the contours of the two-dimensional posterior distribution for Qrms-PS and 
Uq, derived by marginalizing the four-dimensional distribution over Os/i^ and to. These should be 
compared to the corresponding contours for just the smaller-angular-scale data which are shown in 
Fig. 2. The addition of the DMR data significantly tightens the constraints. 

FigTire 7 shows the one-dimensional posterior distribution functions for Qq, derived by marginal- 
izing the four-dimensional ones over the other three parameters. For the flat-A model, using the 
galactic frame DMR data with quadrupole included and corrected for faint high-latitude galactic 
emission, in conjunction with the smaller-angular-scale data, results in — 0.9 being favored 
(first column of Fig. 7), while the ecliptic frame DMR data with quadrupole excluded and no 
further correction for foreground contamination results in = 1 being favored (third column of 
Fig. 7). In the open case, the galactic- frame quadrupole- included DMR data, in conjunction with 
the smaller-angular-scale data, favor f^o ~ 0.5 — 0.6 (second column of Fig. 7), while the ecliptic- 
frame quadrupole-excluded DMR data set combination favor ~ 0.4 (fourth column of Fig. 7). 
Perhaps the most reliable constraints follow from panels k and I of Fig. 7, which are derived from 
the ecliptic frame DMR data with quadrupole excluded and no further correction for foreground 
emission, used in conjunction with the smaller-angular-scale data from the SP94Ka, Python I-III, 
MAX 4 ID, MAX 5 HR, and OVRO experiments. In this case, the favored value of Uo is 1.0 (0.43) 



^*We have checked in a couple of cases that the projected hkeHhood function is fairly similar to the marginalized 
one. This indicates that our "prior" choice of setting to zero the likelihood function outside the parameter range we 
consider does not significantly affect our results. 
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with Q, > 0.17 (0.29 < flo < 0.93) required at 2 a for the flat-A (open) models.^^ Since the two 
= 1 models at the right hand edges of each pair of plots are identical, when the flat-A and open 
model posterior probability distribution functions are renormalized to make this ^Iq = 1 model have 
the same probability in both cases, the probability for = 0.3 (which is close to what is observed, 
see, e.g., Peebles Sz Ratra 2002) is about 50% higher in the open case than in the flat-A case. While 
not highly statistically significant it is still interesting that the open "foil" provides a better fit to 
this data than docs the "favored" flat-A model. It should also be noted that this conclusion also 
follows from the smaller-angular-scale data alone: see the first two columns of Fig. 3. It is puzzling 
why these results are not more consistent with those from less exact analysis of more recent CMB 
anisotropy data, which find that the low-density fiat-A model is favored over the low-density open 
case (see, e.g., Netterfield et al. 2002; Pryke et al. 2002; Stompor et al. 2001). Perhaps the data 
used in this paper has undetected systematic errors, or perhaps it is the more recent data which 
has these (we note that some of the more recent data analyses require use of "prior" assumptions 
before they hone in on the "favored" fiat-A model). It appears that an analysis of more data than 
considered here will be needed to resolve this issue. 

Figure 8 shows the one-dimensional marginalized posterior distribution functions for Qsh"^- 
For both the fiat-A and open models, and for all combined data sets considered, a low il^/i^ ~ 
0.005 — 0.009 is favored. For the most reliable data combination, the DMR ecliptic frame data 
with quadrupole excluded and no other correction for foreground emission, in conjunction with the 
SP94Ka, Python I-III, MAX 4 ID, MAX 5 HR, and OVRO data sets (panels k and I of Fig. 8), 
n^h"^ = 0.007 is favored in both the flat-A and open cases, with Qsh'^ < 0.026 {nsh'^ < 0.027) 
required at 2 cr for the fiat-A (open) models. These values are somewhat lower than what are 
favored by less exact analyses of more recent CMB anisotropy data (see, e.g., Netterfield et al. 
2002; Pryke et al. 2002; Stompor et al. 2001) and those derived from standard nucleosynthesis 
theory and the observed deuterium abundances (Buries, NoUett, Sz Turner 2001); they are, however, 
more consistent with the lower Jlg/i^ value from standard nucleosynthesis theory and the observed 
helium and lithium abundances (Cyburt, Fields, & Olive 2001). These values are consistent with 
those derived using only the smaller-angular-scale data (see the central two columns of Fig. 3) . 

Figure 9 shows the one-dimensional marginalized posterior distribution functions for Iq. The 
flat-A model favors to ~ 12 — 15 Gyr, while the open model favors to ~ 12 — 14 Gyr. The open 
case is very consistent with the results from the smaller-angular-scale anisotropy data alone (see 
the last column of Fig. 3), however, the smaller-angular-scale data alone favors to ~ 18 — 20 Gyr for 
the flat-A model (see the second last column of Fig. 3) , larger than what is favored when the DMR 
data is included in the mix. For the most reliable data combination, the DMR ecliptic frame data 
with quadrupole excluded and no other correction for foreground emission, in conjunction with the 



^^Clearly, the data considered here, which, aside from OVRO, are most sensitive at multipole index I < 200 (OVRO 
is most sensitive at I ~ 600), favor a less steeply rising (with I) spectrum at I < 200, and also favor some significant 
amount of power at / ~ 600, thus favoring an open model over Oo = 1 and the fio = 1 model over a lower density 
flat-A one. 
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SP94Ka, Python I-III, MAX 4 ID, MAX 5 HR, and OVRO data (panels k and I of Fig. 9), to = 13 

Gyr is favored, with to < 19 Gyr required at 2 a for both the open and flat-A cases. These values 
are in excellent accord with estimates for the age of the universe from the ages of the oldest stars 
(see, e.g., Carretta et al. 2000; Krauss &: Chaboyer 2001). 

Figure 10 shows the one-dimensional posterior distribution functions for Qrms-PSj derived by 
marginalizing the four-dimensional ones over the other three parameters. The CMB anisotropy 
data lead to tight constraints on Qrms-PS- 

The peak values of the one-dimensional posterior distributions shown in Figures 7-10 are listed 
in the figure captions for the case when the four-dimensional posterior distributions are normalized 
such that L(Qrins-PS = A^K) = 1. With this normalization, marginalizing over the remaining 
parameter the fully marginalized posterior distributions are, for the fiat-A (open) models: 3 x 
10^^^ (7 X 10^^^) for the galactic-frame quadrupole- included faint-high-latitude-foreground-corrected 
DMR data in conjunction with the smaller-angular-scalc data; 1 x 10*^^^ (5 x 10*^^^) for this data 
less the SuZIE data; 6 x 10^65(7 x lO^^S) for this DMR data and the SP94Ka, Python I-III, MAX 4 
ID, MAX 5 HR, and OVRO smaller-angular-scale data; 3 x 10^" (1 x 10^^) for the ecliptic-frame 
quadrupole-excluded, and no other correction for foreground emission, DMR data in conjunction 
with the smaller-angular-scale data; 2 x 10^^° (1 x 10^-*^) for this data less the SuZIE data; and, 
1 X 10^7^(2 X 10^71) for this DMR data and the SP94Ka, Python I-III, MAX 4 ID, MAX 5 HR, 
and OVRO smallcr-angular-scale data. These numerical values are consistent with the indications 
from Fig. 7 that the open case is mildly favored over the flat-A model. 



7. Conclusion 

We have derived constraints on cosmological model parameters in the open and flat-A CDM 
models, from joint analyses of combinations of the DMR, SP94, Python I-III, ARGO, MAX 4-1-5, 
White Dish, OVRO, and SuZIE GMB anisotropy data sets. The constraints derived here are not 
of very high statistical significance. 

The data considered here mildly favor the open case over the flat-A model, although they do 
not rule out the currently favored flat-A model with ^Iq ~ 0.3. The favored value of in the open 
(J7o ~ 0.4 — 0.6) and flat-A (f^o ~ 0-9 ~ 1) rnodels are somewhat dependent on whether the DMR 
quadrupole moment is included in or excluded from the analysis. Resolving this issue will likely 
require analyses that use new higher quality large angular scale CMB anisotropy and foreground 
emission data. 

Constraints on O^sh^ and to arc only weakly dependent on the data set combination considered. 
The data considered here favors lower ^Ish^ ~ 0.005 — 0.009, or younger, to ~ 12 — 15 Gyr, universes. 

Some of the constraints derived here are mildly inconsistent with those derived elsewhere, 
but gratifyingly not so at high statistical significance. Tighter and more robust constraints on 



- 12 - 



cosmological parameters will require a models-based joint analysis of a larger collection of CMB 
anisotropy data sets. 
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Fig. 1. — Confidence contours and maxima of the two-dimensional posterior probability density 
distribution functions, as a function of the two parameters on the axes of each panel (derived by 
marginalizing the four-dimensional posterior distributions over the other two parameters) . Dashed 
lines (crosses) show the contours (maxima) of the open case and solid lines (solid circles) show 
those of the flat-A model. Contours of 0.25, 0.5, 1, 2, and 3 a confidence are shown (some contours 
are not labelled). Panels a) — c) in the left column show the (rj^/i^, Q,q) plane, while panels d) — /) 
in the right column show the (to, Qq) plane. Panels a) &: d) in the top row are from an analysis of 
all the small scale data. Panels b) & e) in the middle row are from an analysis of all but the SuZIE 
small scale data. Panels c) &; /) in the bottom row are from an analysis of the SP94Ka, Python 
I-III, MAX 4 ID, MAX 5 HR, and OVRO data sets. 
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Fig. 2. — Confidence contours and maxima of the marginalized two-dimensional (Qrms-PS) ^^o) 
posterior probability density distribution functions. Panels a), c), & e) in the left column show 
the flat-A model and panels b), d), & /) in the right column show the open model. Lines show 
the two-dimensional posterior probability density distribution function 1, 2, and 3 a confidence 
limits derived from the various combination smaller-angular-scale data sets considered. Solid circles 
show the maxima of the two-dimensional posterior distributions. Shaded regions show the two- 
dimensional posterior probability density distribution function confidence limits for the DMR data 
(G98; Stompor 1997); densest shading shows the 1 a confidence region and Icss-dcnse shading 
shows the 2 a region. The DMR results are a composite of those from analyses of the two extreme 
data sets: i) galactic frame with quadrupole included and correcting for faint high-latitude galactic 
emission; and ii) ecliptic frame with quadrupole excluded and no other galactic emission correction 
(G98). Panels a) &; b) in the top row are from an analysis of all the small scale data. Panels c) & 
d) in the middle row are from an analysis of all but the SuZIE small scale data. Panels e) &; /) in 
the bottom row are from an analysis of the SP94Ka, Python I-III, MAX 4 ID, MAX 5 HR, and 
OVRO data sets. 
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Fig. 3. — One-dimensional posterior probability density distribution functions for (two left 
most columns), for VtBh^-, (two central columns), and for to (two right most columns), derived by 
marginalizing the four-dimensional ones over the other three parameters, in the open and flat-A 
models. These have been renormalized to unity at the peaks. Dotted lines show the formal 1, 2, 
and 3 a confidence limits derived from these one-dimensional posterior distributions. Panels a), c), 
&L e) in the first column, g), i), &i k) in the third column, and m), o), &l q) in the fifth column show 
the flat-A model and panels b), d), Sz f) in the second column, h), j), k, I) in the fourth column, 
and n), p), and r) in the sixth column show the open model. Panels in the top row are from an 
analysis of all the small scale data, those in the middle row are from an analysis of all but the SuZIE 
small scale data, and those in the bottom row are from an analysis of the SP94Ka, Python I-III, 
MAX 4 ID, MAX 5 HR, and OVRO data sets. When the four-dimensional posterior distributions 
are normalized such that L((5nns-PS = jiK) = 1, the peak values of the one-dimensional 
distributions shown in panels a) - /) are 6 x lO^^S 2 x lO^^^, 3 x lO^'^S 2 x lO^^^, 1 x 10^°^ and 
2 X 10^°^ respectively; in panels g) - 1) are 3 x lO^^^^ 7 ^ iq273^ ^ ^ iq27z ^ g ^ io273^ 5 x 10^03^ and 
7 X 10203, respectively; and in panels m)-r) are 6 x 10^'^°, 1 x 102^^ 3 x lO^^o, 1 x 10^^^ 9 x lO^oo, 
and 1 X 10^°^, respectively. 
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Fig. 4. — One-dimensional posterior probability density distribution functions for Qrms-PS) derived 
by marginalizing the four-dimensional ones over the other three parameters, in the open and fiat-A 
models. Conventions are as described in the caption of Fig. 3. Solid vertical lines show the ±1, 
±2, and ±3 o confidence limits derived by projecting the corresponding small-scale data four- 
dimensional posterior distributions. Also shown are the 2 a DMR (marginalized and projected) 
confidence limits; these are a composite of those from the two extreme DMR data sets (see caption of 
Fig. 2). When the four-dimensional posterior distributions are normalized such that L((5rms-PS = 
/nK) = 1, the peak values of the one-dimensional distributions shown in panels a) — /) are 
9 X 10^'^°, 1 X \{F^, 4 X 10^'^°, 8 X 10^™, 1 X lO^o^, and 1 x lO^o^, respectively. 
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Fig. 5. — Confidence contours and maxima of the two-dimensional posterior probability density 
distribution functions, as a function of the two parameters on the axes of each panel (derived by 
marginalizing the four-dimensional posterior distributions over the other two parameters). Dashed 
lines (crosses) show the contours (maxima) of the open case and solid lines (solid circles) show 
those of the flat-A model. Contours of 0.25, 0.5, 1, 2, and 3 a confidence are shown (some contours 
are not labelled). Panels a) — c) in the first column and panels g) — i) in the third column show 
the (Jlfi/i^, Oo) plane, while panels d) — f) in the second column and panels j) — I) in the fourth 
show the {to, ^o) plane. Panels a) — /) in the two left hand columns use the DMR galactic frame 
data with quadrupole included and corrected for faint high-latitude galactic emission while panels 
g) — I) in the two right hand columns use the DMR ecliptic frame data with quadrupole excluded 
and no other galactic emission correction (G98) . Panels in the top row are from an analysis of the 
DMR and all smaller-angular-scale data, those in the middle row are from an analysis of all but 
the SuZIE data, and those in the bottom row are from an analysis of the DMR, SP94Ka, Python 
I-III, MAX 4 ID, MAX 5 HR, and OVRO data sets. 
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Fig. 6. — Confidence contours and maxima of the marginalized two-dimensional (Qrms-PS) ^o) 
posterior probability density distribution functions. Panels a), c), &: e) in the first column and 
g), i), &: k) in the third column show the flat-A model while panels b), d), & /) in the second 
column and h), j), & /) in the fourth column show the open model. Shaded regions show the 
two-dimensional posterior probability density distribution function confidence limits for the DMR 
data (see Fig. 2 caption for details). Lines show the two-dimensional posterior probability density 
distribution function 1,2, and 3 a confidence limits derived from the various combination data sets 
considered. Solid circles show the maxima of the two-dimensional posterior distributions. Lines 
in panels a) — /) in the two left hand columns use the DMR galactic frame data with quadrupole 
included and corrected for faint high-latitude galactic emission while lines in panels g) — I) in the 
two right hand columns use the DMR ecliptic frame data with quadrupole excluded and no other 
galactic emission correction (G98) . Panels in the top row are from an analysis of the DMR and all 
smaller-angular-scale data, those in the middle row are from an analysis of all but the SuZIE data, 
and those in the bottom row are from an analysis of the DMR, SP94Ka, Python I-III, MAX 4 ID, 
MAX 5 HR, and OVRO data sets. 
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Fig. 7. — One-dimensional posterior probability density distribution functions for Oq, derived by 
marginalizing the four-dimensional ones over the other three parameters, in the open and flat-A 
models. These have been renormalized to unity at the peaks. Dotted lines show the formal 1, 2, 
and 3 a confidence limits derived from these onc-dimensional posterior distributions. Panels a), c), 
& e) in the first column, and g), i), & k) in the third column show the flat-A model while panels 6), 
d), & /) in the second column, and h), j), &l I) in the fourth column show the open model. Panels 
a) — f) in the two left hand columns use the DMR galactic frame data with quadrupole included and 
corrected for faint high-latitude galactic emission while lines in panels g) — I) in the two right hand 
columns use the DMR ecliptic frame data with quadrupole excluded and no other galactic emission 
correction (G98). Panels in the top row are from an analysis of the DMR and all smaller-angular- 
scale data, those in the middle row are from an analysis of all but the SuZlE data, and those in 
the bottom row are from an analysis of the DMR, SP94Ka, Python 1-111, MAX 4 ID, MAX 5 HR, 
and OVRO data sets. When the four-dimensional posterior distributions are normalized such that 
-^(<5rms-PS = A*K) = 1, the peak values of the one-dimensional distributions shown in panels 
a) — I) are 5 x 

6 X 10^'^^ 2 X 10371, and 4 x 103^1, respectively. 
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Fig. 8. — One-dimensional posterior probability density distribution functions for Q.Bh?i derived 
by marginalizing the four-dimensional ones over the other three parameters, in the open and flat-A 
models. Conventions are as described in the caption of Fig. 7. When the four-dimensional posterior 
distributions are normalized such that L(Qrms-PS = = 1, the peak values of the one- 

dimensional distributions shown in panels a) — I) are 2 x 10^^^, 6 x 10^^''', 1 x 10^^^, 3 x lO'^'^''', 
4 X 10^6^, 4 X 10^^^, 3 X 10^''2^ I X 10443^ 2 x lO''^^ 9 x lO^''^^ g ^ io372^ and 1 x 10^^^, respectively. 
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Fig. 9. — One-dimensional posterior probability density distribution functions for to, derived by 
marginalizing the four-dimensional ones over the other three parameters, in the open and flat-A 
models. Conventions are as described in the caption of Fig. 7. When the four-dimensional posterior 
distributions are normalized such that i^(Qrms-PS = = 1, the peak values of the one- 

dimensional distributions shown in panels a) — I) are 4 x 10^^^, 1 x 10^^^, 2 x 10"^^"^, 8 x 10^^^, 
9 X 10^64^ 9 X lO^^"^, 6 X lO''^^, 3 x 10^"^°, 3 x lO'^^^ 2 x lO""^", 2 x 10^™, and 3 x lO^^o, respectively. 
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Fig. 10. — One-dimcnsional posterior probability density distribution functions for Qrms-PS; derived 
by marginalizing the four-dimensional ones over the other three parameters, in the open and flat-A 
models. Conventions are as described in the caption of Fig. 7. Also shown here as solid vertical lines 
are ±1, ±2, and ±3 a confidence limits derived by projecting the corresponding four-dimensional 
posterior distributions, as well as the 2 a DMR (marginalized and projected) confidence limits 
which are a composite of those from the two extreme DMR data sets (see caption of Fig. 2) . When 
the four-dimensional posterior distributions are normalized such that L(Qrms-PS = yuK) = 1, 
the peak values of the one-dimensional distributions shown in panels a) — I) are 
4 X 10-^3^, 1 X 10^^35^ 2 X 10365, 2 X 10^65, 1 X 10^^°, 2 X 10"^"^°, 5 x lO^^g^ 2 x 10'^'^°, 4 x lO^^o, and 
3 x 10^'''°, respectively. 



